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Abstract 

In this article, we show how to cahbrate the widely-used SVI parameterization 
of the implied volatility smile in such a way as to guarantee the absence of static 
arbitrage. In particular, we exhibit a large class of arbitrage-free SVI volatility 
surfaces with a simple closed-form representation. We demonstrate the high quality 
of typical SVI fits with a numerical example using recent SPX options data. 

1 Introduction 

The stochastic volatility inspired or SVI parameterization of the imphed volatihty smile 
was originally devised at Merrill Lynch in 1999 and subsequently publicly disseminated 
in [T^ . This parameterization has two key properties that have led to its popularity with 
practitioners: 

• For a fixed time to expiry t, the implied Black-Scholes variance o"|g(fc, t) is linear in 
the log-strike as \k\ — > oo consistent with Roger Lee's moment formula [23] • 

• It is relatively easy to fit listed option prices whilst ensuring no calendar spread 
arbitrage. 

The consistency of the SVI parameterization with arbitrage bounds for extreme strikes 
has also led to its use as an extrapolation formula [20j. Moreover, as shown in [15], the 
SVI parameterization is not arbitrary in the sense that the large-maturity limit of the 
Heston implied volatility smile is exactly SVI. However it is well-known that SVI smiles 
may be arbitrageable. Previous work has shown how to calibrate SVI to given implied 
volatility data (for example [2Z])- Other recent work [6] has been concerned with showing 
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how to parameterize the volatihty surface in such a way as to preclude dynamic arbi- 
trage. There has been some work on arbitrage-free interpolation of implied volatilities or 
equivalently of option prices [T], [TT], [IB], [21] ■ Prior work has not successfully attempted 
to eliminate static arbitrage and indeed, efforts to find simple closed-form arbitrage-free 
parameterizations of the implied volatility surface are still widely considered to be futile. 

In this article, we exhibit a large class of SVl volatility surfaces with a simple closed- 
form representation, for which absence of static arbitrage is guaranteed. Absence of 
static arbitrage — as defined by Cox and Hobson [8] — corresponds to the existence of a 
non-negative martingale on a filtered probability space such that European call option 
prices can be written as the expectation, under the risk-neutral measure, of their final 
payoffs. This definition also implies (see [H]) that the corresponding total variance must 
be an increasing function of the maturity (absence of calendar spread arbitrage). Us- 
ing some mathematics from the Renaissance, we show how to eliminate any calendar 
spread arbitrage resulting from a given set of SVI parameters. We also present a set of 
necessary conditions for the corresponding density to be non-negative (absence of butter- 
fly arbitrage), which corresponds — from the definition of static arbitrage — to call prices 
being decreasing and convex functions of the strike. We go on to use the existence of 
such arbitrage-free surfaces to devise a new algorithm for eliminating butterfly arbitrage 
should it occur. With both types of arbitrage eliminated, we achieve a volatility surface 
that typically calibrates well to given implied volatility data and is guaranteed free of 
static arbitrage. 



In Section 2.1, we present a necessary and sufficient condition for the absence of 



calendar spread arbitrage. In Section |2.2[ we present a necessary and sufficient condition 
for the absence of butterfly arbitrage, or negative densities. In Section [3} we present 
various equivalent forms of the SVI parameterization. In Section |4| we exhibit a large and 
useful class of SVI volatility surfaces that are guaranteed to be free of static arbitrage. In 
Section [5| we show how to calibrate SVI to observed option prices, avoiding both butterfly 
and calendar spread arbitrages. We further show how to interpolate and extrapolate in 
such a way as to guarantee the absence of static arbitrage. Finally, in Section |6| we 
summarize and conclude. 

Notations. In the foregoing, we consider a stock price process {St)^yQ with natural 
filtration {J-'t)t>oy we define the forward price process {Ft)^yQ by Ft := K{St\J-'o). 
For any A; G M and t > 0, Cbs(^,o'^^) denotes the Black-Scholes price of a European 
Call option on S with strike Fte'', maturity t and volatility a > 0. We shall denote the 
Black-Scholes implied volatility by aBs{k,t), and define the total implied variance by 

The implied variance v shall be equivalently defined as v{k,t) = cr'^s{k,t) = w{k,t)/t. We 
shall refer to the two-dimensional map {k,t) i— w{k,t) as the volatility surface, and for 
any fixed maturity t > 0, the function k h-> w{k, t) will represent a slice. We propose below 
three different — yet equivalent — slice parameterizations of the total implied variance, and 
give the exact correspondence between them. For a given maturity slice, we shall use the 
notation w{k; x) where x represents a set of parameters, and drop the t-dependence. 
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2 Characterisation of static arbitrage 



In this section we provide model-independent definitions of (static) arbitrage and some 
preliminary results. We define static arbitrage for a given volatility surface in the following 
way, which is equivalent to the definition of static arbitrage for call options recalled in the 
introduction (see also [25]). 

Definition 2.1. A volatility surface is free of static arbitrage if and only if the following 
conditions are satisfied: 

(i) it is free of calendar spread arbitrage; 
(a) each time slice is free of butterfly arbitrage. 

In particular, absence of butterfly arbitrage ensures the existence of a (non-negative) 
probability density, and absence of calendar spread arbitrage implies monotonicity of 
option prices with respect to the maturity. The following two subsections analyse in 
details each of these two types of arbitrage, in a model-independent way. 



2.1 Calendar spread arbitrage 



Calendar spread arbitrage is usually expressed as the monotonicity of European call option 
prices with respect to the maturity (see for example [3] or [S]). Since our main focus 
here is on the implied volatility, we translate this definition into a property of the implied 
volatility. Indeed, assuming proportional dividends, we establish a necessary and sufficient 
condition for an implied volatility parameterization to be free of calendar spread arbitrage. 
This can also be found in [llj and [13] and we outline its proof for completeness. 

Lemma 2.1. // dividends are proportional to the stock price, the volatility surface w is 
free of calendar spread arbitrage if and only if 

dtw{k, t) > 0, for allk eR and t > 0. 



Proof. Let (Xt)^>Q be a martingale, L > and < ti < t2- Then the inequality 

E[(X,,-L)+] >E[(Xi, -L)+] 

is standard. For any i = 1, 2, let Cj be options with strikes Ki and expirations tj. Suppose 
that the two options have the same moneyness, i.e. 

Then, if dividends are proportional, the process {Xt)t>Q defined by Xt := St/ Ft for all 
t > is a martingale and 



9i 



{Xt^-e'^) 



k\ + 
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So, if dividends are proportional, keeping the moneyness constant, option prices are non- 
decreasing in time to expiration. The Black-Scholes formula for the non-discounted value 
of an option may be expressed in the form CBs{k,w{k,t)) with Cbs strictly increasing 
in its second argument. It follows that for fixed k, the function w{k,-) must be non- 
decreasing. □ 



Lemma 2.1 motivates the following definition. 

Definition 2.2. A volatility surface w is free of calendar spread arbitrage if 

dtw{k, t) > 0, for allk eR and t > 0. 



2.2 Butterfly arbitrage 



In Section 2.1 , we provided conditions under which a volatility surface could be guaranteed 



to be free of calendar spread arbitrage. We now consider a different type of arbitrage. 



namely butterfly arbitrage (Definition 2.3). Absence of this arbitrage corresponds to the 
existence of a risk-neutral martingale measure and the classical definition of no static 
arbitrage, as developed in [12j or In this section, we consider only one slice of the 
implied volatility surface, i.e. the map k h-> w{k,t) for a given fixed maturity t > 0. For 
clarity we therefore drop — in this section only — the t-dependence of the smile and use 
the notation w{k) instead. Unless otherwise stated, we shall always assume that the map 
k w{k, t) is at least of class C^(]R) for all t > 0. 

Definition 2.3. A slice is said to he free of butterfly arbitrage if the corresponding density 
is non-negative. 

Recall the Black-Scholes formula for a European call option price: 

CBs{k,w{k)) = S {Afid+{k)) - e^J\f{d^{k))) , for all A; e M, 



where Af is the Gaussian cdf and d±{k) := —k/ ^Jw{k) ± y/w{k)/2. Let us define the 
function (7 : M — M by 



,„):.A_^^V_^f 1 H.1V^. (2.1) 



2w{k) ) 4 V^(^) ' 4 

This function will be the main ingredient in the determination of butterfiy arbitrage as 
stated in the following lemma. 

Lemma 2.2. A slice is free of butterfiy arbitrage if and only if g{k) > for all k E M. 
and lim d+{k) = —00. 

A;— >+oo 
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Proof. It is well known [2] that the probability density function p may be computed from 
the call price function C as 



p{k) 



d^c{k) 



d^CBs{k,w{k)) 



for any k G 



K=Fte^ 



K=Fte'= 

Explicit differentiation of the Black-Scholes formula then gives for any k G 

gik) f d.{kr 



p{k) 



^2TTw{k) 



exp 



We have so far implicitly assumed that the function p is a well-defined density, and in 
particular that it integrates to one. This may not always be the case though, and one 
needs to impose asymptotic boundary conditions. In particular, call prices must converge 
to as /c tends to infinity, which is equivalent to having limfc_>_|_oo dj^{k) = — oo. We refer 
the reader to [21] for a proof of this equivalence. □ 



3 SVI formulations 

We first recall here the original SVI formulation proposed in [13], and then present some 
alternative (but equivalent) ones. We emphasize in particular that even though the origi- 
nal ( "raw" ) formulation is very tractable and has become popular with practitioners, it is 
difficult — seemingly impossible — to find precise conditions on the parameters to prevent 
arbitrage. 



3.1 The raw SVI parameterization 

For a given parameter set xr = {(^yb, p,m,a}, the raw SVI parameterization of total 
implied variance reads: 

u^ik; Xr) = a + b |p (A; - m) + \/ {k - m)^ + a^j , (3.1) 

where a G M, 6 > 0, |p| < 1, m G M, a > 0, and the obvious condition a + b a a/1 — > 0, 
which ensures that w{k\XR) ^ for all A; G M. This condition indeed ensures that 
the minimum of the function w{-;xr) is non-negative. Note further that the function 
k I—)- w{k\ Xr) is (strictly) convex on the whole real line. It follows immediately that 
changes in the parameters have the following effects: 

• Increasing a increases the general level of variance, a vertical translation of the smile; 

• Increasing b increases the slopes of both the put and call wings, tightening the smile; 

• Increasing p decreases (increases) the slope of the left (right) wing, a counter-clockwise 
rotation of the smile; 
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• Increasing m translates the smile to the right; 

• Increasing a reduces the at-the-money (ATM) curvature of the smile. 

We exclude the trivial cases p — I and p = —1, where the volatility smile is respectively 
strictly increasing and decreasing. We also exclude the case a — which corresponds to 
a linear smile. 



3.2 The natural SVI parameterization 

For a given parameter set xn — 1^, P,^X}j ^^e natural SVI parameterization of total 
implied variance reads: 

w{k; ;^^) = A + 1 1 1 + Cp (A; - //) + ^(C(^ - //) + + (1 - , (3.2) 

where a; > 0, A G M, yU G M, |p| < 1 and ( > 0. It is straightforward to derive the 
following correspondence between the raw and natural SVI parameters: 

Lemma 3.1. We have the following mapping of parameters between the raw and the 
natural SVI: 



{a, b, p, m,a)^ |^A + | (l - p^) ^ ^ _ -J—JLj ^ (3.3) 

and its inverse transformation, between the natural and the raw SVI: 

(A...P...C) . („ - I (1 - /) -^.p. -Jt^. . (,4) 

3.3 The SVI Jump- Wings (SVI-JW) parameterization 

Neither the raw SVI nor the natural SVI parameterizations are intuitive to traders in 
the sense that a trader cannot be expected to carry around the typical value of these 
parameters in his head. Moreover, there is no reason to expect these parameters to 
be particularly stable. The SVI- Jump- Wings (SVI-JW) parameterization of the implied 
variance v (rather than the implied total variance w) was inspired by a similar parame- 
terization attributed to Tim Klassen, then at Goldman Sachs. For a given time to expiry 
t > and a parameter set xj — {''^t, i^tiPti Ct, ^t} the SVI-JW parameters are defined from 
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the raw SVI parameters as follows: 



Vt 



+ b |— pm + y/m'^ + a^} 



, t 

Vt =^fe(l-p), (3-5) 

1 



^* = ^ (^a + & o- V 1 - p2 

with := fft. Note that this parameterization has an explicit dependence on the time 
to expiration t, and hence can be viewed as generalizing the raw (expiration-independent) 
SVI parameterization. The SVI-JW parameters have the following interpretations: 



• Vt gives the ATM variance; 

• gives the ATM skew; 

• 'Pt gives the slope of the left (put) wing; 

• Q gives the slope of the right (call) wing; 

• Vt is the minimum implied variance. 



If smiles scaled perfectly as Ij^/wl (as is approximately the case empirically), these 
parameters would be constant, independent of the slice t. This makes it easy to extrapolate 
the SVI surface to expirations beyond the longest expiration in the data set. Also note 
that by definition, for any t > we have 

Wt = 



dk 



k=0 



The choice of volatility skew as the skew measure rather than variance skew for exam- 
ple, reflects the empirical observation that volatility is roughly lognormally distributed. 
Specifically, following the lines of p3l Chapter 7] , assume that the instantaneous variance 
process satisfies the SDE 

dvt = a{vt) dt + rjy/ut f3{vt) dZt, for alH > 

where t] > 0, {Zt)t>o is a standard Brownian motion and a and /3 two functions on M+ 
ensuring the existence of a unique strong solution to the SDE (see for instance [2^ for 
exact conditions), then the ATM variance skew 



t^o dk 



k=0 
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exists and is proportional to /3(f). If the variance process is lognormal so that /5(f) 
behaves hke \/f , the hmit of the at-the-money volatility skew as time to expiry tends 
to zero is constant and independent of the volatihty level. This consistency of the SVI- 
JW parameterization with empirical volatility dynamics thus leads in practice to greater 
parameter stability over time. The following lemma provides the inverse representation 



of (3.5). 



Lemma 3.2. Assume that m 7^ 0. For any t > 0, define the (t- dependent) quantities: 



13 ■■= P 1 and a := sign(/3) J— - 1. 



where we have further assumed that /3 G [—1, iQ Then, the raw SVI and SVI- J W param- 
eters are related as follows: 



P = 1 



(ct + Pt) 

VtxfWt 



a = Vft — ba^yl — p"^, 

{Vt - Vt) t 



m 



b |— p + sign(Q;) vT+~Q^ — «a/1 — 
a m. 

If m = 0, then the formulae above for b, p and a still hold, but a = {vtt — a) /b. 



a = am. 



Proof. The expressions for 6, p and a follow directly from (3.5). Assume that m 7^ and 



let j3 := p — 2ijjt\/uh/b and a := a/m G M. Then the expressions in (3.5) give 

sign (a) 



f3 



which implies that 



a = sign(/3)J^-l. 



Using (3.5), we also have 

(ft - Vt) t 
b 

from which we deduce m in terms of a, and the expression of cr is recovered from the 



|— p + sign(a) vT+~Q^ — a a/1 — p^ | 



equality a = am. The expression for a in the case m = is straightforward from (3.5 ). □ 



^The condition /? e [—1, 1] is equivalent to —pt < 2-04 < ct, i.e. to the convexity of the smile. 
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3.4 Arbitrage and absence thereof in SVI parameterizations 



Given a volatility surface, it is natural to wonder whether it is free of arbitrage. Since 
we can easily switch from any of the SVI formulations to either of the other two using 
Lemma 3.3| and Lemma 3^, we shall state the following results only for the raw SVI pa- 
rameterization (3.1). Referring to (3.1) as a volatility surface is a slight abuse of language 



since (3.1) is really an expiry-independent slice parameterization. A volatility surface is 



thus understood as a (discrete) collection of slices, with a different set of parameters for 



each expiry. Checking calendar arbitrage in the sense of Lemma 2.1 is then equivalent to 



checking for calendar arbitrage for any pair of expiries ti and t2- The following lemma 
establishes a sufficient condition for the absence of calendar spread arbitrage. 



quartic polynomial (given in (3.7) below) has no real root. 



Lemma 3.3. The raw SVI surface (3.1) is free of calendar spread arbitrage if a certain 



Proof. By definition, there is no calendar arbitrage if for any two dates ti ^ t2, the corre- 
sponding slices w (■, ti) and w (■, ^2) do not intersect. Let these two slices be characterised 
by the sets of parameters xi '■= {^i, &i, pi, mi} and X2 '■= {^2, &2, c"2, P2, '"^2}, and as- 
sume for convenience that < ti < ^2- We therefore need to determine the (real) roots 
of the equation w {k,ti) = w {k,t2)- The latter is equivalent to 



ai + biS pi {k - mi) + y {k - mi) + (tI> = 02 + 62 S P2 (^ - "^2) + 



{k - 7712) + a2 

(3'.6) 



Leaving \/{k — mi) + af on one side, squaring the equality and rearranging it leads to 
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(« + /3fc) ^J{k-m2)^ + al = bj{ik- mi)' + a^} -bl{{k- m2f + al] - (« + I3ky 



where a := 02 — ai + bipimi — 62P2"^2 and (3 := 62P2 — &iPi- Squaring the last equation 
above gives a quartic polynomial equation of the form 



k'^ + as k^ + a2 k'^ + ai k + ao = 0, 



(3.7) 



where each of the coefficients lengthy yet explicit expression^ in terms of the parameters 
{oi, 61, pi, CTi, mi} and {02, &2, P2, c"2, ■'^2}- If this quartic polynomial has no real root, then 
the slices do not intersect and the lemma follows. Roots of a quartic polynomial are known 
in closed-form thanks to Ferrari and Cardano [3] . Thus there exist closed-form expressions 
in terms of Xi and X2 for the possible intersection points of the two SVI slices. □ 



Remark 3.1. // the quartic polynomial (3.7) has one or more real roots, we need to 



check whether the latter are indeed solutions of the original problem (3.6), or spurious 



solutions arising from the two squaring operations. The absence of real roots of the quartic 
polynomial is clearly a sufficient — but not necessary — condition. 



^Explicit expressions for these coefficients can be found in the R-code posted on http: //faculty. 
|baruch . cuny . edu/jgatheral^ 
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Remark 3.2. By a careful study of the minima and the shapes of the two slices w{-,ti) 
and w{-,t2), it is possible to determine a set of conditions on the parameters ensuring no 
calendar spread arbitrage. However these conditions involve tedious combinations of the 
parameters and will hence not match the computational simplicity of the lemma. 

For a given slice, we now wish to determine conditions on the parameters of the raw 
SVI formulation (3.1) such that butterfly arbitrage is excluded. By Lemma 



2.1 



this is 

equivalent to showing (i) that the function g defined in (2.1) is always positive and (ii) 



that call prices converge to zero as the strike tends to infinity. Sadly, the highly non-linear 
behaviour of g makes it seemingly impossible to find general conditions on the parameters 
that would eliminate butterfly arbitrage. We provide below an example where butterfly 
arbitrage is violated. Notwithstanding our inability to find general conditions on the 
parameters that would preclude arbitrage, in Section |4} we will introduce a new sub-class 
of SVI volatility surface for which the absence of butterfly arbitrage is guaranteed for all 
expiries. 



Example 3.1. (From Axel Vogt on wilmott.com) Consider the raw SVI parameters: 



{a,b,m,p,a) = (-0.0410,0.1331,0.3586,0.3060,0.4153) 



(3.^ 



with t = 1. These parameters give rise to the total variance smile w and the function g 
(defined in (2.1)) on Figure^ where the negative density is clearly visible. 




-1.5 -1.0 -0.5 0.0 0.5 




Figure 1: Plots of the total variance smile w (left) and the function g defined in (2.1) 



(right), using the parameters (3.8) 



4 Surface SVI: A surface free of static arbitrage 



We now introduce a class of SVI volatility surfaces — which we shall call SSVI (for 'Surface 



SVr) — as an extension of the natural parameterization (3.2). For any maturity t > 0, 
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define the at-the-money (ATM) implied total variance 6t '■= (T^^{0,t)t. We shall assume 
that the function 6 is at least of class on M^. An ATM option with zero time to expiry 
has no value so 6*0 := limt_^o 6'* = 0. 

Definition 4.1. Let ip be a smooth function from to M.*^ such that the limit limj_j.o 6t(p{6t] 
exists in M. We refer to as SSVI the surface defined by 

w{k,et) = ^^\i + pm)k + ^ {^{et)k + + {1 - p^) I . (4.1) 



2 

From Section |3| SSVI corresponds to the natural SVI volatility surface parameter- 



ization (3.2) with xn = {0, 0, p, ^t, (/9(6't)}. Note that this representation amounts to 
considering the volatility surface in terms of ATM variance time, instead of standard 
calendar time, similar in spirit to the stochastic subordination of [7]. 



Remark 4.1. In the parameterization (4.1), the ATM variance curve 9t may be viewed 



as a (vector) parameter of the volatility surface. Moreover, this parameter is directly 
observable given market prices for a finite set of expiries, and can be considered well- 
known to traders even for expiries which are not explicitly quoted. The explicit reference 
to 6t also emphasizes the importance of studies such as J 10] of the ATM variance structure 
in classical models which may shed some light on how to impose dynamics on SSVI. 

The ATM implied total variance is 6t = (j|s(0,t)t and the ATM volatility skew is 
given by 

5feCTBs(fc,t)|fc=0 = 7As^9kW{k,et 

Furthermore the smile is symmetric around at-the-money if and only if p = 0. This 
is consistent with [H Theorem 3.4] which states that in a standard stochastic volatility 
model, the smile is symmetric if and only if the correlation between the stock price and 
its instantaneous volatility is null. Since = 0, we have at time t = 0: 

w{k,eo) = ]^<i)Q {pk+\k\) , for any A; gM, (4.3) 

where 0o '■= hnie-s-o 0(p{9). 0o = is characteristic of stochastic volatility models as in Ex- 



= ^^'W). (4.2) 



ample 4A 0o > as in Example |4 . 2 1 gi ves a V-shaped time zero smile which is characteris- 
tic of models with jumps and in particular, characteristic of empirically observed volatility 
surfaces. For notational convenience, we shall always assume that YrnvtyAoodt = oo. As 
proved in [2^, this is equivalent (assuming no interest rate) to the stock price (assumed 
to be a non-negative martingale) to converging to zero as t tends to infinity. Although 
this holds in many popular models (Black-Scholes, Heston, exponential Levy), this is not 
always true, see [12] for counter-examples. If lim^ 9t is finite, all our results remain 
valid, but only on the support of the function t Of. 

The following theorem gives precise necessary and sufficient conditions to ensure that 



the SSVI volatility surface (4.1) is free of calendar spread arbitrage (Lemma 2.1) and 
also matches the term structure of ATM volatility and the term structure of the ATM 
volatility skew. 
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Theorem 4.1. The SSVI surface (4.1) is free of calendar spread arbitrage if and only if 



1- dtOt > 0, for all t > 0; 

2. < deiOifie)) <j,{l + yr^) ^{0), for all 6 > 0, 
where the upper bound is infinite when p = 0. 



In particular, this theorem imphes that the SSVI surface (4.1 ) is free of calendar spread 
arbitrage if the skew in total variance terms is monotonically increasing in trading time 
and the skew in implied variance terms is monotonically decreasing in trading time. In 
practice, any reasonable skew term structure that a trader defines has these properties. 

Proof. Since the definition of calendar spread arbitrage does not depend on the log- 
moneyness k, there is no loss of generality in assuming k fixed. First note that dtw{k, 9t) = 



dgw{k,6t)dt6t so the SSVI volatility surface (4.1) is free of calendar spread arbitrage if 
dew{k, e)>0 for all ^ > 0. 

Consider first the case |p| < 1. To proceed, we compute, for any 9 > 0, 

2dew{k, 9) = ipo{x, p) + 7(6')V^i(x, p), 
with X := k^{9), 7(^) := de{9i^{9)) /ip{9), 

, , \ 1+ px \ 1 X + p 

ip(}{x,p) := 1 H 1^^^=^^^= and ipi{x,p) := x 



^/x^T2pxTl ' ' [ y/x^ + 2pX + 1 

For any |p| < 1, ?/'o(a;,p) is strictly positive for all x G M. Now define the set 

( (-00,0) U (-2p,oo), ifp<0, 
Vp=l (-oo,-2p) U (0,oo), ifp>0, 

[r\{o}, ifp = 0. 

Then ^i(-, p) > if x € I^p and V^i(-, p) < if x G M \ {V^ U {0, -2p}). It follows that 

Mx,P) 



dewik, 9)>0ii and only if { J^nt' ^ (4.4) 

7W<-7V^' for xGM\(I?pU {0,-2 p}). 



When X G {0, — 2p}, then ijji{x, p) = and so dew{k,9) > 0. The inequalities (4.4) thus 
give necessary and sufficient conditions for absence of calendar spread arbitrage for any 
given X G M. To determine the tightest possible bounds on 7(6*), we compute 

V^o(a;,p) n ^ -f V^o(a;,p) 1 + ^/l - p^ 

sup = and ml — 



xe©p V'i(a;,p) xGR\(CpU{o-2p}) V^i(x,p) p2 
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The supremum in the first equahty is never attained (the function increases to zero from 
below as tends to infinity). However the infimum in the second equality is attained at 
X = —p ^ Vp. It follows that 



dew{k, ^) > if and only if < -f{9) < 



1 + yT^ 



Note that when p = 0, the infimum above is taken over an empty set, and there is hence 
no upper bound. 



When p = 1, for any [k, 6*) G M x (0, oo), we have 

1 + x 



dew{k, 6) 



1 + 



X] 



[l + l{0)x) 



2(l + 7(^)a;) if x > -1, 
otherwise. 



Obviously, dgw{k,9) > if a; > 0. For x > —1, clearly dgw{k,9) > if and only if 
7(0) G [0, 1]. Similarly, with p = —1, we have 



dow{k, 9) 



1 — X 



[l-^x) 



2{l--f{9)x) ifa;<l, 
otherwise. 



Again dgw{k, ^) > if a; < 0, and for x < 1, dgw{k, ^) > if and only if -f{9) G [0, 1]. □ 



The following lemma is a straightforward consequence of (3.3) and (3.5). 



Lemma 4.1. The S VI- J W parameters associated with the SS VI surface (4.1) are 



Pt 

Ct 



Ot/t, 

\pV0t^{9t) 
1 



2 
1 

2 ' 
9t 



/9tv{9t){l-p), 
/9,^{9,){l + p), 



We now give several examples of SSVI implied volatility surfaces (4.1). 

Example 4.1. A Heston-like parameterization 

Consider the function Lp defined by 



\9 



1 -e- 



with A > 0. Then for all 9 > 0, we immediately obtain 



A202 



-xe 



+ \9-l 
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For any A > 0, the map 6 ^ de {9ip{9)) /'{>{6) is strictly decreasing on (0, oo) with limit as 

6 tends to zero equal to one. Since the quantity (1 + ■>/! — p^)/ is greater than one for 
any p e [—1, 1], the conditions of Theorem 4-l\ o-re satisfied. This function is consistent 
with the implied variance skew in the Heston model as shown in fT^ Equation 3.19]. 

Example 4.2. Power-law parameterization 

Consider ^{6) = r^Q-^ with r] > and < < 1. Then de {O^iO)) /(piO) = 1 - 7 e (0, 1) 
holds for all 6 > 0, and hence the conditions of Theorem o^^e satisfied. In particular if 

7 = 1/2 then Lemma 4-i implies that the SVI- J W parameters ipt, Pt, o-nd Ct associated with 



the SSVI volatility surface (4.1) are constant and independent of the time to expiration t. 



Furthermore, Equation \4.2\ implies that the ATM volatility skew is given by 



dkCrBs{k,t)\ 



k=0 



PV 
2Vi' 



The following theorem provides sufficient conditions for a SSVI surface (4.1 ) to be free 
of butterfly arbitrage. 



Theorem 4.2. The SSVI volatility surface (4.1) is free of butterfly arbitrage if the fol- 
lowing conditions are satisfied for all 9 > 0: 

1. M^)(1 + IPI)<4; 

2. ^^(^^)2(l + |p|) <4. 

Proof. For ease of notation, we suppress the explicit dependence of 9 and (p on t. By 
symmetry, it is enough to prove the theorem for < p < 1. We shall therefore assume 



so, and we define z := (pk. The function g defined in (2.1) reads 



64 (^2 + 2zp + 1) 



3/2- 



where 



f{z) := a - hv'9 - -y9\ 
lb 



and where a, h and c depend on z. In the following, we frequently use the inequality 

z^ + 2zp + I = [z + pf + I - p^ > 0. 
Computing the coefficient of ip'^9'^ in f{z) explicitly gives 

c= ^Jz^ + 2zp+l\^[l + p") (z + p)' + 2p(z + p)v/^2^2^pTI+ (l-p') p'} 
> ^z^ + 2zp+l { (1 + p2) {z + pf + 2p{z + pf + (1 - p2) p2} 
= + 2zp + 1 {(1 + pf {z + pf + (1 - p2) p2} >0. 
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Thus if 



we have 



0<Oip< 



4 4 
and < ^v?^ < 



1 + P 
46 



a — 



/(^)> 



(1+P)' 

It is then straightforward to verify that 
2/i(^) 



1 + p' 

if6>0, 
if 6 < 0. 



(1 + P)^ 



= + 2;zp + 1 {z^p - 2(1 - p)p + 2(1 + p) (1 - p') + p} 

+ p (^ + p) V 3p (1 - p2) + 2 (1 - p2) - ^p (z^ + 2zp + 1) , 



which is clearly positive for z < 0. To see that fi{z) is also positive when z > 0, we 
rewrite it as 

2/i(^) 



(1 + p)^ 

z2 + 22p + l-(z + p)}|p(;.-i^y + 2(l + p) (l-p')+p(l-^^^ 

+ (l + p){^(2-p2)+2(l + p) {l-p^)+p}. 
Consider now the function f2{z). It is straightforward to verify that 

f2iz) = --0^ + + 2^P + 1)'^' + 2{z^ + 2zp + 1) + ^z^ + 2zp + l 

which is positive by inspection if 2; < 0. To see that f2{z) is also positive when 2; > 0, we 
rewrite it as 



f2iz) = Z 



3 1 + P^ 



+ 3z^p + 2{z^ + 2zp + 1) 



(1 + P)' 

+ {z^ + 2zp + l) ^^^/z^ + 2zp+l-{z + p)^ 
+ ^Jz'^ + 2zp + l + 2 V + Z + p. 



Thus f{z) > in all cases. From Lemma 2.2 , we are left to prove that limfc_>oo d^{k) 



00. 



A straightforward computation shows that this is satisfied as soon as Condition 1 in The- 
orem 14.21 holds. □ 



Remark 4.2. A SSVI volatility surface (4.1) is free of butterfly arbitrage if 



-Jn^ max (pt, q) < 2, and {pt + Ct) max (pt, q) < 2, 
hold for all t > 0. The proof follows from Lemma 4J_ by re-expressing Conditions^ and 



^ of Theorem 4-2 in terms of S VI-, JW parameters. 
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The following lemma shows that Theorem |4.2| is almost if-and-only-if . 



Lemma 4.2. The SSVI volatility surface (4.1) is free of butterfly arbitrage only if 

6^(6) (1 + IpI) < 4, for all 6 > 0. 
Moreover if 6ip{6) (1 + |p|) = 4, the SSVI surface is free of butterfly arbitrage only if 

M^)'(i + |p|)<4. 

Thus Condition [I] of Theorem 



4-2 is necessary and Condition is tight. 



Proof. Considering the SSVI surface (4.1) and the function g defined in (2.1), we have 
16-^VW(1 + P)' 



4-Mg)Mi + p) 

64 ^ 8f{e)k 

16-^VW(1-P)' 9^(9)^1 -p) 



64 



^^{9)k 



+ 
+ 



as A; — 7- +00, 
as A; —J- — oo. 



The result follows by inspection. 



□ 



Remark 4.3. The asymptotic behavior of SSVI (4.1) as \k\ tends to infinity is 

(1 ± P) 9t 



w{k,9t 



-ip{9t)\k\ + 0{l), foranyt>0. 



We thus observe that the condition 9(p{9) (1 + |p|) < 4 o/ Theorem 4-2 corresponds to the 
upper bound of 2 on the asymptotic slope established by Lee 123] and so again, Condition^ 
of Theorem \4.2\ is necessary. 



The following corollary follows directly from Theorems |4. 1| and 4.2 



Corollary 4.1. The SSVI surface (4.1) is free of static arbitrage if the following condi- 
tions are satisfied: 

1- dt9t > 0, for allt>0 

2. < dei9ip{9)) <j2{l + \A^) ip{9), for all 9 > 0; 

3. 9ip{9) (1 + IpI) < 4, for all 9 > 0; 

4. 9ip{9y (1 + IpI) < A, for all 9 > 0. 



Remark 4.4. Consider the function (^{9) = rj9^^ with i] > from Example 4-2, then 
Condition^ imposes 7 G (0,1). From Condition^ such surfaces can be free of static 
arbitrage only up to some maximum expiry. Take for instance the simple case 9t := cr'^t 
for some a > 0. Then the map ip : t ^ 9tip{9t) (1 + |p|) — 4 is clearly strictly increasing 
with ip{0) = —4 and limt_5.oo ip{t) = 00. Therefore there exists tQ>0 such that ipit) < 
for t < t*Q. The map ■ t ^ 9t'{){9tf (1 + |p|) - 4 zs 
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strictly increasing if 'j & (0, 1/2) with V'2(0) = —4 and lim ilj{t) = +oo; t/iere exists 
tl>0 such that ip2it) <0 fort<tl. 



• strictly decreasing if j E (1/2,1) with lim?/'2(0) = +oo and lim ip{t) = —4; there 
exzsfo > such that ■il)2it) < for t>t\. 

• constant if a = 1/2 with %l)2 = —4. 

When'-) G (0, 1/2), i/ie surface is guaranteed to he free of static arbitrage only fort < tQ/\t\. 
For 7 G (1/2, 1), this remains true only for t G (0,^9) fl {t\,oo) (which may he empty). 
When 7 = 1/2, static arbitrage cannot occur fort < t^. However, the behavior for large 9 
can be easily modified so as to ensure that the entire surface is free of static arbitrage. For 
example, the choice 

gives a surface that is completely free of static arbitrage provided that (1 + |p|) < 2. 



Remark 4.5. In the Heston-like parameterization of Example note that 

lim M^^)(l + |p|) = ^-^. 
Therefore Condition^ of Corollary 4-1 imposes A > (1 + |p|) /4. 



The following model-independent theorem provides a way to expand the class of volatil- 
ity surfaces that are guaranteed to be free of static arbitrage by adding a suitable time- 
dependent function. 



Theorem 4.3. Let {k,t) 1— )■ w{k,t) be a SSVI volatility surface (4.1) satisfying the con 



ditions of Corollary 4-1 (in particular free of static arbitrage), and a : IR+ — ?■ IR+ a non- 
negative and increasing function of time. Then the volatility surface {k,t) 1— )■ Wa{k,6t) '■= 
w{k, 6t) + at is also free of static arbitrage. 



Proof. Since dtWa{k, 6t) := dtw{k, 6t) + dtat, Lemma [23] implies that Wa is free of calendar 
spread arbitrage if dtat > and at > 0. We now show that Wa is also free of butterfly 
arbitrage. For clarity, since butterfly arbitrage does not depend on the time parameter t, 
we shall use the simplified notation w{k) := w{k,9t), and likewise Waik) := Wa{k,9t). 



Similarly, in view of (2.1), we shall define the map ga{k), where the function w is replaced 
by Wa- We consider the case p < since the case p > follows by symmetry, and the 
result is obvious when p = 0. Let us consider the function Ga : M — ?► M defined by 

Ga{k) := g{k) - ga{k), for all k eR, 

and let k* := —2p/ip{9t) > be the unique solution to the equation w'{k) = 0. We can 
compute explicitly the following: 

Ga{k) = ^ (^r- - -^Tt] f 4fc + W'{k) - W'{k)e (^r- + -^Tt] ] , 

4 \Waik) w{k)J\ \Waik) w{k)JJ' 
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which imphes 



w'jk) {w'{k)+ Ak) Waik) - 2k^w'{k) 
4 



(4.6) 



Wa{ky 

Since w'{0) = p9tf{9t) < the equation w'{k) + 4A; = has a unique solution k^: > 0, and 
w'{k) + Ak is strictly positive for any k > k^, and strictly negative when k < k^. By strict 
convexity of the function w it also follows that k^ < k*. Therefore for any k G {k^,,k*), 
the two inequalities w'{k) < and w'{k) + 4:k > hold, and therefore d^Gaik) > 0. Since 
by construction Go{k) = 0, we therefore conclude that g{k) > ga{k) for any k e (/c*, k*). 
For k ^ (fc*, k*), the inequality ^^(/c) < ga{k) holds as soon as daGa{k) < 0. Consider first 



the case k > k*. We can rewrite (4.6) as 



daGa{k) 



w'{k) 2k [2w^{k) - kw'{k)] + wMw'{k) 
4 WaikY 



so that it suffices to prove the inequality 2wa{k) — kw'{k) > for any k > k*. It suffices 
to prove daGa{k) < for then we have the inequality gaik) > g{k) > and there is no 
butterfly arbitrage. 

First consider the case k > k*, so that w'{k) > 0. Recall that a continuously differen- 
tiable function / is convex on the interval (a, b) if and only if f{x) — f{y) > f'{x){x — y) 
for all (x, G (a, 6). Setting x = k and y = 0, we conclude that 2wa{k) — kw'{k) > 
since Wa(0) > 0. It follows that daGa{k) < for any k > k*. 

For any A; < 0, we always have w'{k) < 0, the inequality 2wa{k) — kw'{k) > follows 
by convexity as above, and hence daGa{k) < for any A; < 0. We prove here that 
Qaik) > 5'a(0) for all such k. Since we already showed that (7a(0) > 0, the result follows. 



From the definition of ga and (|2.1|) 

9a{k) - gaiO) = 



kw'{k) 



2{w{k) + a) 
w'{kf [ 1 



4 



w{k) + a 
w"{0) 



- 1 
1 

^4 



w{0) + a 



(4.7) 



A straightforward analysis shows that the function k H- w"{k) is strictly increasing on the 
interval {0,k*/2) and strictly decreasing on {k*/2,k*). The easy computation w"{0) = 
w"{k*) implies that w"{k) > w"{0) on (0, k*). Also, w'{OY > w'{kf on (0, k*). Simplifying 



(|4/rj), it follows that 



kw'{k) 



> 



1 

4 



2{w{k) 



a 



1 + 



w'{ky 



w{0) + a w{k) + a J w{k) + a 



w{0) + a 
kw'{k) 



w'jkf 
w{k) + a 
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Note that w'{kf < w'{0) w\k) < w'{Of on the interval (0, k*) so 



gaik) - gaiO) > 



w'{0)w'{k) 



k w'{k) 



(4i 



4 V'^(O) + ^ w{k) + a J w{k) + a 
We now prove the following claim: kw{<d) - - w(0)] > for A; G (0, A;*). Indeed, 



1 - 



Condition 2 of Theorem 



4.2 



implies that 1 



ek 



^Jip'^k'^ + 2^pk + l. 



> 0. Then (recall that p < 0) the 



right-hand side of the above equality represents an increasing function on (0, k*) which is 
equal to zero at the origin, and the claim holds. Then, from ( 4.8[ ), 



-w'{k) 



(w(0) + a) {w{k) + a) 



k {w{0) + a) - [w{k) - w{0)] 



> 0. 



□ 



Remark 4.6. Given a set of expirations < ti < . . . < tn (n > 1) and at-the-money 
implied total variances Q < 9t^ < . . . < 9t^, Corollary 4-i gives us the freedom to match 
three features of one smile (level, skew, and curvature say) but only two features of all the 
other smiles (level and skew say), subject of course to the given smiles being themselves 
arbitrage- free. Theorem \4.3\ may allow us to match an additional feature of each smile 
through at- 



5 Numerics and calibration methodology 
5.1 How to eliminate butterfly arbitrage 

In Section |4| we showed how to define a volatility smile that is free of butterfly arbitrage. 
This smile is completely defined given three observables. The ATM volatility and ATM 
skew are obvious choices for two of them. The most obvious choice for the third observable 
in equity markets would be the asymptotic slope for k negative and in FX markets and 
interest rate markets, perhaps the ATM curvature of the smile might be more appropriate. 



In view of Lemma 4.1[ supposing we choose to fix the SVI-JW parameters Vt, ipt and pt 



of a given SVI smile, we may guarantee a smile with no butterfly arbitrage by choosing 
the remaining parameters and v[ as 

Ct = + 2 i)u and v^ = Vt -3 . 

In other words, given a smile defined in terms of its SVI-JW parameters, we are guaranteed 
to be able to eliminate butterfly arbitrage by changing the call wing Cf and the minimum 
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variance Vt, both parameters that are hard to cahbrate with available quotes in equity 
options markets. 



Example 5.1. Consider again the arbitrageable smile from Example \3.1\ The correspond- 
ing S VI- J W parameters read 

(vt, tpt, Pt, ct, Vt) = (0.01742625, -0.1752111, 0.6997381, 1.316798, 0.0116249) . 

We know then that choosing {ct,Vt) = {c°,v°) := (0.3493158,0.01548182) gives a smile 
free of butterfly arbitrage. It follows by continuity of the parameterization in all of its 
parameters, that there must exist some pair of parameters {dl,v^) with G (c°,q) and 
vl G {ytiV") such that the new smile is free of butterfly arbitrage and is as close as 
possible to the original one in some sense. In this particular case, choosing the objective 
function as the sum of squared option price differences plus a large penalty for butterfly 
arbitrage, we arrive at the following '^optimal" choices of the call wing and minimum 
variance parameters that still ensure no butterfly arbitrage: 

{c;,v;) = (0.8564763,0.0116249). 

Note that the optimizer has left Vt unchanged but has decreased the call wing. The resulting 
smiles and plots of the function g are shown in Figure 111 





Figure 2: Plots of the total variance smile (left) and the function g defined in (2.1 ) (right) 



using the parameters (3.8). The graphs corresponding to the original Vogt parameters is 



solid, to the guaranteed butterfly-arbitrage-free parameters dashed, and to the "optimal" 
choice of parameters dotted. 



Remark 5.1. The additional flexibility potentially afforded to us through the parameter 
at of Theorem sadly does not help us with the Vogt smile of Example \5.1 . For at 
to help, we must have at > 0; it is straightforward to verify that this translates to the 
condition Vt (1 — p^) < Vt which is violated in the Vogt case. 
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5.2 Calibration of SVI parameters to implied volatility data 



There are many possible ways of defining an objective function, the minimization of which 
would permit us to calibrate SVI to observed implied volatilities. Whichever calibration 
strategy we choose, we need an efficient fitting algorithm and a good choice of initial 
guess. The approach we will present here involves taking a square-root fit as the initial 
guess. We then fit SVI slice-by-slice with a heavy penalty for calendar spread arbitrage 
(i.e. crossed lines on a total variance plot). Consider two SVI slices with parameters Xi 
and X2 where ^2 > ^i- We first compute the points ki {i = 1, . . . ,n) with n < 4 at which 
the slices cross, sorting them in increasing order. If n > 0, we define the points ki as 

ki := ki — 1, 

ki := ^{ki^i + ki), if 2 < i < n. 
For each of the n + 1 points ki, we compute the amounts q by which the slices cross: 



max 



0,w{ki;xi) -w{ki;x2] 



Definition 5.1. The crossedness of two SVI slices is defined as the maximum of the Ci 
{i = 1, . . . ,n) . If n = 0, the crossedness is null. 

An example SVI calibration recipe 

• Given mid imphed volatilities aij = crBs(fcj, tj), compute mid option prices using 
the Black-Scholes formula. 

• Fit the square-root SVI surface by minimizing sum of squared distances between 
the fitted prices and the mid option prices. This is now the initial guess. 

• Starting with the square-root SVI initial guess, change SVI parameters slice-by 
slice so as to minimize the sum of squared distances between the fitted prices 
and the mid option prices with a big penalty for crossing either the previous 



slice or the next slice (as quantified by the crossedness from Definition 5.1) 



There are obviously many possible variations on this recipe. The objective function 
may be changed and when finally working to optimize the fit slice-by-slice, one can work 
from the shortest expiration to the longest expiration or in the reverse order. In practice, 
working forward or in reverse seems to make little difference. Changing the objective 
function on the other hand will make some difference especially for very short expirations. 

5.3 Interpolation and extrapolation of calibrated slices 

Suppose we follow the above recipe above to fit SVI to options with a discrete set of 
expiries. In particular, each of the resulting SVI smiles will be free of butterfly arbitrage. 
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It's not immediately obvious that we can interpolate these smiles in such a way as to 
ensure the absence of static arbitrage in the interpolated surface. The following lemma 
shows that it is possible to achieve this. 

Lemma 5.1. Given two volatility smiles w{k,ti) and w{k,t2) with ti < t2 where the two 
smiles are free of butterfly arbitrage and such that w{k, T2) > w{k, Ti) for all k, there exists 
an interpolation such that the interpolated volatility surface is free of static arbitrage for 
ti <t <t2. 



Proof. Given the two smiles w{k,ti) and w{k,t2), we may compute the (undiscounted) 
prices C{Fi, Ki,ti) =: Cj of European calls with expirations {i = 1,2) using the Black- 
Scholes formula. In particular, since the two smiles are free of butterfly arbitrage, 

^>0, for. = 1,2. 

Consider any monotonic interpolation 6t of the at-the-money implied total variance w{0, t). 
Let Ki = FiQ^ and Kt = Fte^. Then for any ti < t < t2, define the price Ct = C{Ft, Kt, t) 
of a European call option to be 

— = at— + {l-at)—, (5.1) 
Kt Ki K2 

where for any t G (^1,^2), we define 

V - V 



By construction, for fixed k, the inequality 

>0 

holds so that there is no calendar spread arbitrage. Also, because of the square-roots in 



drKt 



the definition (5.2), the at-the-money interpolated option price will be almost perfectly 
consistent with the chosen implied total variance interpolation 6t. Moreover, if the two 
smiles w{k, ti) and w{k, t2) are free of butterfly arbitrage, we have dK,KC{k, t) > 0. To see 



this, first note that because all the options have the same moneyness, the identity (5.1) 
is equivalent to 

-^ = at-^ + {l-at)-^. (5.3) 
rt ri r2 

Then note that the ratio C{F, K,t)/F is a function of F and K only through the log- 
moneyness k. Also, for K = Kt, Ki, K2, we have 

dk^ dk ' 
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Applying this to (5.3), we obtain 



All the terms on the rhs are non-negative, so the Ihs must also be non-negative. We 
conclude that there is no butterfly arbitrage in the interpolated slice and thus that there 
is no static arbitrage. The interpolated volatility surface may be retrieved by inversion of 
the Black-Scholes formula. □ 



We could conceive of a myriad of algorithms for extrapolating the volatility surface. 
For example, one way to extrapolate a given set of n > 1 (arbitrage-free) volatility smiles 
with expirations < ti < . . . < t„ would be as follows: At time to = 0, the value of a call 
option is just the intrinsic value. We may then interpolate between and ti using the 



algorithm presented in the proof of Lemma 5.1, thereby guaranteeing no static arbitrage. 



For extrapolation beyond the final slice, we suggest to first recalibrate the final slice using 



the SSVI form (4.1 ). Then fix a monotonic increasing extrapolation of 6t (asymptotically 
linear in time would seem to be reasonable) and extrapolate the smile for t > tn according 
to 

wik,et) = wik,etj + et-et„, 



which is free of static arbitrage if w{k,6t^) is free of butterfly arbitrage by Theorem 4.3 



5.4 A calibration example 

We take SPX option quotes as of 3pm on 15-Sep-2011 (the day before triple- witching) and 
compute implied volatilities for all 14 expirations. The result of fitting square-root SVI is 
shown in Figure |3j The result of fitting SVI following the recipe provided in Section 5.2 



is shown in Figure |4j With the sole exception of the first expiration (options expiring at 
the market open on the following morning), the fit quality is almost perfect. 
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Figure 3: Red dots are bid implied volatilities; blue dots are offered implied volatilities; 
the orange solid line is the square-root SVI fit 

6 Summary and conclusion 

We have found and described a large class of arbitrage-free SVI volatility surfaces with a 
simple closed-form representation. Taking advantage of the existence of such surfaces, we 
showed how to eliminate both calendar spread and butterfly arbitrages when calibrating 
SVI to implied volatility data. We have also demonstrated the high quality of typical SVI 
fits with a numerical example using recent SPX options data. The potential applications 
of this work to modelling the dynamics of the implied volatility surface are left for future 
research. 
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Figure 4: Red dots are bid implied volatilities; blue dots are offered implied volatilities; 



the orange solid line is the SVI fit following recipe of Section 5.2 
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